A possible chiral spin-liquid phase in non-centrosymmetric i?BaCo407 
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Based on a symmetry approach, we propose a possible explanation of the weak ferromagnetic 
component recently observed in YBaCosFeOr (Valldor et al. Phys Rev B, 84 224426 (2011)) and 
other isostructural compounds in the high-temperature spin-liquid phase. Due to the polar nature 
of their crystal structure, a coupling between time-odd scalar spin chirality which we suggest as 
the primary order parameter and macroscopic magnetization is possible as follows from the general 
form of the appropriate free-energy invariant. The deduced pseudoproper coupling between both 
physical quantities provides a unique possibility to study the critical behaviour of the chiral order 
parameter. 
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The concept of symmetry and its spontaneous break- 
ing at a phase transition is a central topic of condensed 
matter physics. Invar iance of physical properties of crys- 
tals under spatial operations (such as rotations and re- 
flections) and time-reversal is the basis for their com- 
prehensive description and classiflcation. Time-reversal 
symmetry is the crucial ingredient to understand proper- 
ties and describe phase transitions in magnetic systems 
with long-range spin-dipole ordering. Combination of 
this symmetry with the spatial operations form a pow- 
erful tool of magnetic symmetry based on which such 
fundamental phenomena as antisymmetric exchange 
piezomagnetism [2] and linear magnetoelectric effect 
have been predicted. 

Besides conventional magnets exhibiting spontaneous 
breaking of time-reversal invariance at magnetic phase 
transitions, there is another type of systems where long- 
range spin ordering is suppressed by geometrical frus- 
tration. Usually, the geometrical frustration is caused 
by a combination of triangular-based exchange topol- 
ogy and antiferromagnetic interactions between neigh- 
bour spins. These systems fluctuate between degenerated 
ground states in the so called "spin-liquid" phase, and 
can stay in this regime down to very low temperatures. 
The absence of long-range dipole magnetic ordering does 
not necessarily mean that spin-liquids keep time-reversal 
invariance. Some other spin quantities such as toroidal 
moment (T) and scalar spin chirality (x) have been pro- 
posed as possible candidates for the time-reversal break- 
ing order parameters as well [4-7]. In spite of consider- 
able efforts to flnd systems with actual realization of the 
spontaneous long-range toroidal or chiral order, in the 
absence of dipole ordering, only very recently a possible 
candidate of the latter has been found after observation of 
anomalous Hall effect in the metallic pyrochlore Pr2lr207 

ii. 

In insulating systems, however, identiflcation of the 
long-range chiral ordering is more challenging but is not 



hopeless due to a possible coupling of the chiral order pa- 
rameter to other measurable physical quantities. In this 
letter, we discuss this possibility for the geometrically 
frustrated antiferromagnets i?BaMe407 (i?=rare earth 
or Y, Me=Co,Fe) with a polar hexagonal lattice. Based 
on a simple symmetry consideration, we show that the 
polar nature of the crystal structure allows an appropri- 
ate coupling of the spin-chiral order parameter to sponta- 
neous magnetization (M) along the hexagonal axis. We 
propose this phenomenon as an explanation for the weak 
ferromagnetic component recently observed by Valldor et 
al. [10] in the spin-liquid YBaCosFeOy. 

The prototype crystal structure of the YBaCosFeOy 
composition like in other members of the i?BaMe407 
family has hexagonal Pdsmc symmetry pjj and consists 
of corner-sharing Me04 tetrahedra arranged in alternat- 
ing Kagome and triangular layers stacked along the c 
axis (Fig. [T]) [l2|, [l3|. The Me-sublattice forms a new 
exchange topology where geometrical frustration is in- 
herent to the presence of both trigonal bipyramids and 
triangular clusters [13, EH • 

In spite of extremely strong exchange interactions in 
YBaCosFeOy (paramagnetic Curie- Weiss temperature is 
^ -2000K), no long-range magnetic ordering has been ob- 
served down to IK [10]. It should be pointed out that 
some of the i?BaCo407 compounds exhibit spin ordering 
below ~ lOOK, which has been shown to be closely related 
to a strong magnetoelastic coupling reducing the crystal 
structure symmetry down to monoclinic P2i [l6|. Spin 
dynamics probed by high-resolution inelastic neutron 
scattering revealed a freezing below T/ ~50K [10], thus 
in a wide temperature range above Tj, YBaCosFeOy is 
in a spin-liquid state. Surprisingly, a very weak ferro- 
magnetic component (~ 10~^/iB per Co/Fe ion) along 
the c axis and a hysteretic behaviour of the magnetiza- 
tion in magnetic fleld have been revealed below 590K. 
The existence of the small spontaneous magnetization 
in polycrystalline samples with the YBao.5Cao.5Co407 
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TABLE I. Matrix of irreducible representations (irreps) for 
generators of PGsmcl' space group associated with wave vec- 
tor k=0 and physical quantiles transforming according to 
these representations, t is the time-reversal operator. 
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FIG. 1. (Color online) Schematic representation of the hexag- 
onal PGsmc crystal structure of i?BaMe407 as a corner- 
shared Me04 tetrahedral network. Metals ions in Kagome 
and triangular layers are denoted as Mel (red) and Me2 (yel- 
low), respectively. Solid lines represent exchange bonds. 



(below 387K [I3) and YBaCo407 (at room tempera- 
ture [HI, [lH) compositions had been reported before as 
well, indicating that this remarkable phenomenon is not 
a unique property of the Fe doped crystal. However, the 
higher transition temperature and the fact that the weak 
ferromagnetic component is observed in very clean single 
crystals and is anisotropic makes the YBaCosFeOy case 
a more compelling one. These observations along with 
numerous neutron diffraction measurements proving the 
lack of a coherent magnetic Bragg scattering above lOOK 
113, [HI, EqI, Ell directly indicate that the time-reversal 
symmetry in the pure i?BaCo407 or Fe doped composi- 
tions is broken without dipole spin ordering. The fun- 
damental question immediately arises, what breaks the 
time-reversal invariance in all these materials? 

As it was mentioned above, a spin-dipole order- 
ing is not the only order parameter which breaks time- 
reversal symmetry. Hypothetically, toroidal moment or 
scalar chirality are also considered as the possible vari- 
ants of the time-odd order parameters. The latter is a 
three-spin quantity defined as a mixed product of spins 
X = Si • [Sj X Sk] in a triangular plaquette. The scalar 
chirality takes a nonzero value for non-coplanar spin 
states and its sign represents whether the non-coplanar 
structure is either right- or left-handed. In spite of the 
time-odd nature of both T and x order parameters, they 
have different transformational properties and therefore 
different types of coupling to macroscopic quantities. 
From symmetry point of view, the spin chirality defined 
above is a pseudoscalar i.e. it changes sign at a mir- 
ror reflection and stays invariant under rotations as it is 
demonstrated in Fig. [21 where a result of reflection by 



the mirror planes and the tree-fold rotation of the P6smc 
space group are shown. In group theoretical language, it 
means that the chiral order parameter transforms as the 
time-odd irreducible representation of the PGsmcV 
space group (see Table U we use Miller and Love notation 
for the irreducible representations as it is implemented 
into ISOTROPY program [22|, the superscript indicates 
parity with respect to time inversion). A ferromagnetic 
component along the 6-fold axis (M^) also transforms 
accordingly this representation and can therefore be lin- 
early coupled to X- Taking into account that magnetiza- 
tion is a vector and F^ (g) F^ 3 and polarization along 
the c axis (P^) transforms as F]^, the general form of the 
free-energy invariant relating the macroscopic quantities 
can be written as x^P- Thus, scalar chirality can be 
coupled to magnetization only in polar crystals and in 
a general case, where polarization direction is not fixed 
by symmetry, a condensation of the chiral order param- 
eter must induce magnetization in the direction of the 
polarization. In the case where polarization is frozen like 
in YBaCosFeOr (or other members of the family), cou- 
pling between x and M is pseudoproper and the critical 
behaviour of both order parameters is identical. This 
offers the unique possibility to directly measure the criti- 
cal exponent of the chiral order parameter with sufficient 
accuracy to verify the universality class of this kind of 
phase transitions. 

In-plane (Tx^Ty) and out-of-plane (T^) components of 
the macroscopic toroidal moment transform as the time- 
odd Fg and F^ representations, respectively (Table [Tj) 
and cannot form an appropriate coupling invariant lin- 
ear in respect of M^. 

The absence of both spatial and time inversions cre- 
ates sufficient conditions for a linear magnetoelectric ef- 
fect. The magnetic point group in the chiral phase of 
i?BaCo407 is Qm'm' (magnetic space group is PGsm^c^) 
which allows an = a22 and ass components of the mag- 
netoelectric tensor to be non-zero. Thus, measurement 
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FIG. 2. (Color online) Transformational properties of three 
non-coplanar spins in a triangle (for simplicity, we arbitrary 
took these spins to create equal angles between each other 
but actually the angles can be different). The scalar chi- 
rality, changes sign under the mirror reflection, m, and 
time-inversion, t, but does not change it under the rotation 
C3. Intersections of the mirror planes of the triangle with 
the plane containing it are shown as dot-lines. The signs of 
component for each spin are shown as "+ "or "-" 



of the magnetization in electric field or polarization in 
magnetic field can be an additional probe of a long-range 
chiral order in insulating crystals. 

The presented above symmetry-based phenomenolog- 
ical consideration is not able to reveal the microscopic 
origin of the phase transition to the chiral phase. How- 
ever, let us discuss symmetry aspects of the possible in- 
teractions favouring non-coplanar spin arrangement in 
the case of well defined local correlations. First, we note 
that antisymmetic Dzyaloshinskii-Moriya (DM) interac- 
tions can in some circumstances create conditions for a 
non-zero scalar chirality in a spin triangle as it shown 
in Fig. [3] (left). The phenomenological free-energy term 
responsible for these interactions is written in the form 
of Dij • [Si X Sj] where Dy is the so called Dzyaloshinskii 
vector whose direction depends on the crystal symmetry 
[2^. In particular, in the hexagonal lattice of i?BaCo407, 
it must be parallel to the mirror planes bisecting ij bonds 
and can be decomposed into two components; out-of- 
plane Dout (parallel to the c axis) and in-plane Din (per- 
pendicular to the c axis). Due to the polar nature of the 
structure, the out-of-plane component is always either up 
or down in the whole crystal. This component removes 
the degeneracy between the two planar spin configura- 
tions with different vector chirality (defined as a vector 
product of the two neighboring spins, averaged over three 
spin pairs j) I^i ^ ^j]' summation is taken over 

three pairs of sites along the sides of the triangle in an 




FIG. 3. (Color online) Two variants of coplanar spin con- 
figurations in a triangle with parallel (tt) cind antiparallel 
(ti) directions of the vector chirality and Dout component. 
The latter was arbitrary chosen to be in the up-direction (a). 
Impact of the in-plane Din component on these spin configu- 
rations (b). Spin configuration forbidding (left) and allowing 
(right) z-component and scalar chirality to be non-zero (c). 



ant i- clockwise direction) shown in Fig. [3^. Then, the 
Din component deviates the spins from the {ab) plane 
for the case when [Si x Sj] tt Dout and leaves 

them coplanar when [^i ^ ^j] ti Dout- Thus, de- 

pending on the direction of Dout, the antisymmetric DM 
term can select non-coplanar spin configurations. 

The case, when Din does not couple a z-component 
to the spins in a triangle (Fig. [SJd (right)), is of par- 
ticular interest. It is easy to see that admixture of the 
out-of-plane spin component does not change the sym- 
metry of the planar spin configuration and therefore can 
be added to it (Fig. [3t (right)). The coupling, how- 
ever, cannot be linear since the z-component is invariant 
in respect of the three-fold axis, while the planar spin 
configuration is not. So the out-of-plane component has 
the higher symmetry (3m') and can be coupled only as 
a linear-cubic term between the associated order param- 
eters due to the necessary invariance in respect of the 
time inversion. Thus, interactions which relate the in- 
plane and out-of-plane spin components in this case are 
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TABLE II. Symmetrized combinations of the atomic spins, 
S, on the WyckofF position 6c (Kagome layer) transforming 
as components of the two-dimensional order parameter, 
r?i(i = l,2) H. 
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distinct from DM term coupling orthogonal components 
with identical transformational properties (bilinear cou- 
pling). 

More generally, the global distortions of the type 
shown in Fig. [3^ (right) keeping translational invari- 
ance of the crystal can be presented as combinations of 
basis functions of the F^ irreducible representation of 
the PQ'^mcV space group (Table ITT|). The associated two- 
dimensional order parameter {r]i^r]2) forms a general cou- 
pling invariant with the out-of-plane ferromagnetic com- 
ponent, 7^^772 — Mzr]2i which is reduced down to the 
simple linear-cubic term, M^T^f^ {^•>V2) direction. 

The latter is the equilibrium order parameter for the 
most symmetric variant of the spin structure with the 
vector chirality shown in Fig. [3^ (right). Minimization 
of the free-energy part, F = jM^ + PMzr]2i respect of 
M^, gives the non-zero equilibrium value for the latter, 
Mz = Thus, in spite of the passive role of the DM 

interactions, the macroscopic magnetization and scalar 
chirality are expected to be non-zero due to the existence 
of the linear-cubic type of the coupling. The chiral free- 
energy term, x^P? (which shows that the system gains 
energy in the case of the non-coplanar spin arrangement) 
can clearly provide this coupling with some " adopted" to 
local symmetry meaning of the M and P quantities and 
in a general case it represents antisymmetric interactions 
distinct from DM. These interactions imply a "chiral de- 
pendent" part of the exchange Hamiltonian, tending to 
bring three spins to the mutually orthogonal orientation. 
The local meaning of M and P adopted for a single tri- 
angular unit in a structure should ensure vanishing of 
the coupling invariant for some spin configurations whose 
symmetry forbids non-coplanar states as in the case of 
3m point group shown in Fig. [3b (left). In some partic- 
ular cases, the DM exchange can contribute to the free- 



energy chiral term too (as it was shown above), since the 
existence of the bilinear coupling invariant between order 
parameters means existence of a linear-cubic one as well. 
To further clarify the relevant interactions responsible 
for the chiral state in i?BaMe407, numerical simulations 
based on the Monte- Carlo algorithm with different terms 
in magnetic Hamiltonian including the chiral-dependent 
XMP one are desirable. 

In conclusion, based on symmetry consideration, we ar- 
gue that the weak ferromagnetic component, measured 
in YBaCosFeOy and some other isostructural composi- 
tions, in spin-liquid regime, is related to a long-range 
chiral ordering. The coupling between both macro- 
scopic physical quantities is possible in polar crystals 
due to a phenomenological trilinear invariant of the form 
xMP. This invariant, in a general case represents in- 
teractions distinct from Dzyaloshinskii-Moriya ones and 
provides a linear-cubic coupling between orthogonal spin- 
components. Due to the frozen nature of the polariza- 
tion in the hexagonal lattice of i?BaCo407, the coupling 
between x and M is pseudoproper and the critical be- 
haviour of both order parameters is identical. This offers 
a unique possibility to identify the universality class of 
these exotic phase transitions. The obtained result here 
provides a physical ground for a search of other possible 
systems where spin chirality is the primary order param- 
eter. 
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